Full QED+QCD low-energy constants through reweighting 
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The effect of sea quark electromagnetic charge on meson masses is investigated, and first results 
for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks 
is incorporated in quenched QED+fuU QCD lattice simulations by a reweighting method. The 
reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic 
method on 2 + 1 flavor dynamical domain- wall quark ensembles. 



So far most lattice QCD simulations have been per- 
formed neglecting the electromagnetic (EM) charges. In 
order to calculate physical quantities to high precision, it 
is quite important to include and control this contribu- 
tion. Toward this goal, several attempts regarding this 
issue have been done using quenched QED [lH5[- 

Chiral perturbation theory (ChPT) provides an effec- 
tive guide to extrapolate to the physical quark mass point 
in a lattice calculation, and in combination with QED, is 
a powerful tool for addressing isospin breaking. The QED 
effect can be included in the ChPT framework: partially 
quenched ChPT (PQChPT) with QED was first derived 
by Bijnens and Daniclsson Q in the SU{3) flavor ba- 
sis up to next-to-leading order (NLO) and was recently 
extended by some of us to the SU{2) flavor-|-kaon ba- 
sis [3]. In PQChPT, sea and valence quarks are sepa- 
rately treated. Here, we specify the two valence quarks 
in mesons by indices 1 and 3 and the three sea quarks 
(m, d, s) by indices 4-6, and then introduce quark masses 
rrii and quark EM charges qt in units of the fundamental 
EM charge e. The fundamental charge can be set sepa- 
rately between the sea and valence quark sectors, which 
we specify by eg and ev, respectively. Using this notation, 
NLO SU{S) PQChPT tells us the sea EM charge contri- 
bution to the pseudo-scalar (PS) meson mass-squared is 
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with Xij = Boim, +mj), Qs{3) = diag(g4, 95, gg)- M is an 
energy scale below which the effective theory accurately 
describes the theory (QED-I-QCD), and Bq, Fq, C and Yi 
are low-energy constants (EEC's). Determination of Yi 
requires eg 7^ 0, which is not accessible in quenched QED. 
Note that the EEC's generally depend on masses and 
EM charges of heavier dynamical quarks than u, d and 



s. In the following, three- flavors of dynamical quarks are 
assumed: u, d and s. We also mention that a remarkable 
feature in the three-flavor theory. 
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makes many terms vanish and leads to the simple form 
in Eq. ([T]). Recently, RBC and UKQCD coUaborations 
pointed out that the SU{2) ChPT is preferable to the 
SU{?>) ChPT even in the three-flavor dynamical QCD 
simulation, since the s quark mass is not small enough (3] . 
In this case, the sea EM contribution to a pion mass- 
squared, for example, can be written as 
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where Xi = '^B^rui, Qs{2) = diag((3'4, q^). We remark that 
additional EEC's, F/, F/', J, J', K and K' to Eq. ^ 
arise due to lack of the property ^ in the SU{2) case, 
and the EEC's in Eq. ([3]) generally depend on a mass and 
an EM charge of the s quark. 

While the full QED effect can be incorporated in the 
Monte Carlo (MC) evolution of the gauge field config- 
uration, the usual gauge ensemble has been generated 
only with dynamical QCD. However, the full QED ef- 
fect, in principle, can be included using a reweighting 
method, which was first proposed and tested by Dun- 
can, Eichten and Sedgewick The main purpose of 
this work is to show the practicality of the reweighting 
method for incorporating the sea quark EM charge on a 
domain-wall fcrmion (DWF) ensemble originally gener- 
ated with = and the feasibility of obtaining the full 
QED EEC's. (Some applications of reweighting to a re- 
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alistic QED+QCD simulation were recently reported in 
Refs. 

Full QED+QCD (fQED+fQCD) includes a U{1) pho- 
ton field A in addition to the usual SU (3) link variable 
U for the gluon field and fermion field ip. In order to il- 
lustrate the reweighting method, we consider the systern 
with a fermion action Sf[ip, ip, U] = —ipD[U]ilj, where U 
is the combined SU{3) x U{1) gauge link variable asso- 
ciated with a quark with EM charge ge; 



U = U X 
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Here, we assume the photon fields are generated by a 
non-compact U{1) photon action; 

Suii) [^] = I E E i9,A^(x) ~ d^A,{x)f . (5) 



An expectation value for some observable O in 
fQED-|-fQCD is formally related to the one in quenched 
QED+fuU QCD (qQED+fQCD), in which the photon 
fields in the quark determinants are neglected, via 

ln\ (wO)qQED+fQCD If.. 
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introducing a reweighting factor 
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Because the fine structure constant of QED is small, 
Q!EM = e^/(47r) ^ 1/137, an expansion of w in photon 
fields and small quark masses is instructive. Defining a 
sea quark mass matrix Afs(3) for three-flavor dynamical 
quarks, 

Ms{3) = =aiag{m4,m5,me), m= ,(8) 

m 6 

the reweighting factor is expanded in sea EM charge eg 
and sea quark mass m; 

w[U, U] = l + 'u;itr(A^s(3)Qs(3))™es + W2trQl^^^e'^ 

+0{rn^es,rnele^,), (9) 

with coefficients w^'s, and we used, again, the remark- 
able feature in the three-flavor case ([2]) , making the lead- 
ing contribution in the expansion vanish. Note that 
while the remaining leading contribution in Cg is formally 
0{es), there exists a suppression factor from the small- 
ness of the quark masses and SU{3) flavor symmetry, 
tr(Ms(3)Qs(3))r?T/AQCD, which leads to a comparable or 
smaller contribution than the 0{e^) terms. 

In this study a stochastic estimate with random Gaus- 
sian noise vectors are used to calculate the determinants 
in the reweighting factor. In general, the error due to 
this estimation could be severely underestimated because 
the distribution of the estimator may be asymmetric 
with a long tail. To overcome this difficulty, breaking 
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FIG. 1. Normalized reweighting factor w[U, U] with the EM 
charge Cs = e on each gluon configuration. 



up the determinant into many small pieces is efficient, 
since the effects of the outliers are largely suppressed 
il[i3|. For the splitting, we use a mathematical iden- 
tity for the determinant, so called the n"^-root trick: 
w ~ det r2 = (det fi^/") , which is easily implemented 
by the rational approximation [llj. We apply reweight- 
ing to 2 -}- 1 flavor dynamical DWF and Iwasaki gluon 
configurations generated by the RBC-UKQCD collabo- 
rations The configuration set is one of the ensem- 
bles used in the quenched QED study (3], whose sim- 
ulation parameters are ;5qcd = 2.13, x T x Lg = 
16^x32x16, inverse lattice spacing a^^ = 1.784(44) GeV, 
{am^,amd,amg) = (0.01,0.01,0.04). The U{1) photon 
fields, which have been already generated in the quenched 
QED study, are combined with the gluon configurations 
according to Eq. We also employ n = 24 roots and 
use 4 complex random Gaussian noise vectors per root 
on each configuration to estimate the reweighting factors. 
Fig. [T] shows the obtained factors normalized by the con- 
figuration average. The fiuctuation among configurations 
is moderate, controlled within a factor of 5. 

DWF's explicitly break chiral symmetry due to finite 
size Lg in the extra 5th dimension which can be quan- 
tified by an additive, residual, quark mass for each fla- 
vor. The residual mass of the ensemble in this study 
is a?Tiros(QCD) = 0.003148(46) (in the chiral limit). The 
quenched QED studies Ql show that the valence EM 
charges further shift the quark mass by an amount of 



C('3EMa"T'rcs(QCD))- The Same effect also arises from the 
sea quark charges. This lattice artifact induces a term 
like e2,5rostrQ2^3) in the S'C/(3) ChPT formula ^. (Sim- 
ilar modifications are also needed in the SU (2) formula 
Q.) Here we measure the sea EM charge contribution 
to the residual mass and subtract it from AA/pg. 

In practical simulations where the number of gauge 
configurations is finite, contributions arise from "hair" , 
or photon emission to, and absorption from, the vacuum 
which averages to zero in the large ensemble limit. In 
Ref. it was shown that this hair is a large source 
of noise in hadron correlators. The leading unwanted 
piece can, however, be removed by averaging over plus 
and minus EM charges, the so-called ±e trick, and it 
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FIG. 2. Effective mass for the meson, ami = am3 = 
0.01 (lower) and 0.02 (upper). The blue symbols represent 
the effective masses after applying the reweighting, while 
the red symbols depict the ones without the reweighting. 
The fit results of the masses with uncorrelated fit in 
t are denoted by the horizontal lines. In fitting the full 
QED data, xVd-o.f. [uncorr] = 0.11(Iower)/0.13(upper) and 
X^/d.o.f.[corr] = 0.67(lower)/0.76(upper). 



provides a great advantage in which the unphysical noise 
is exactly canceled in the valence sector 0, |^ ; 



{0{+e,) + 0{-e,)}^Oiel), 
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where 0(ev) represents some observable with a valence 
EM charge Cv There is also "hair" in the sea sector. To 
remove the leading contribution from both the sea and 
valence sectors, we use an averaging, 



1 



{0(+e„ +ev) + 0{-es, -e^)} = 0(e^, e^Cv, e^), (11) 



in the reweighting. Note that the noise from hair associ- 
ated with Cs is already small by virtue of Eq. ([2]) in the 
three-flavor case. 

Using the reweighting factor obtained in this work and 
the meson correlators in the quenched QED study 3, 
the reweighted meson correlators are obtained by Eq. ^ . 
Examples of effective masses for the 7r+ meson are shown 
in Fig. [2] For the fit results of the masses, we take the 
same fit range (t ~ 9 — 16) as in Ref. Q and also perform 
both correlated (corr) and uncorrelated (uncorr) fits 
in t. Although the reweighting tends to push up the 
tt'^ masses, the gap is quite hard to detect outside of 
statistical errors. 

With the reweighted data of the meson masses calcu- 
lated, chiral fits are performed to obtain the QED LEC's 
in Eqs. ([!]) and ([3|). Although C is known from the 
quenched QED study [1], it provides a valuable consis- 
tency check with the quenched QED result. In fitting 
for the LEC's, we anticipated a problematic hierarchy 
between the and egBv terms, attributable to a double 



TABLE I. Transformation property under Eqs. (fT3)) - (fT5)) for 
individual terms in NLO 5(7(3) and 5(7(2) PQChPT. 



transformation 


terms in NLO PQChPT associated with 
Yi, Yl, Y{' C, J, J' K., IC' 


Ti (Eq. (Ha}) 
T2 (Eq. dM}) 
% (Eq. din}) 


even even even 
even odd odd 
even even odd 
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which left the egGv terms unresolved. Although the dif- 
ficulty can, in principle, be overcome with enormous 
statistics, drastic improvements are provided by engi- 
neering sign flips in the EM charge. Besides the ±e trick 
(Eqs. ([TU| and (fTTj) '). consider a basic transformation 



Ti : (mi,(7i;m3,g3) — > (ma, 173; mi, gi), 

under which the meson system is invariant (CPT 
addition to Ti, let us introduce transformations: 



T2 ■■ (mi,(7i;m3,g3) 
73 : (mi,(?i;m3,g3) 



(mi, 

(^3, 



-qi; mi. 



-93) 
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Eqs. ([T3|) - ()15p form a set of transformations that ex- 
change two valence quark masses and EM charges with, 
or without, flipping the sign of ev Note that T2 and Ts 
yield only partial invariances of Eqs. dl]) and dS]), in the 
sense that the invariance holds only for specific terms in 
each. In Tab.|ll the transformation property of each term 
in NLO PQChPT is summarized. While the and egCv 
terms retain their even and oddness under 7i and T2 to 
all orders in quark mass, the transformation property un- 
der Ts is not preserved at order higher than 0{am) in the 
quark mass expansion. At NLO in SU{2) PQChPT in 
formula dSJ, the esCv term is a mixture of even and odd 
contributions since the three-flavor feature d21) is explic- 
itly broken. By adding and subtracting squared meson 
masses related by these transformations, each term can 
be separately extracted. 

Figs. [3]|8] show individual sea-quark charge contribu- 
tions to the pion mass-squared, Cg, esev(73-even) and 
esev(73-odd) parts. The lattice artifact ingredient, which 
is caused by the finiteness of L^, is subtracted from the 
el term. In the figures, we can clearly see that the hi- 
erarchy between the el and egCv term is 0(10^), as ex- 
pected by the suppression given by Eq. (|12[) , and the sep- 
aration using the transformation T2 successfully works. 
The valence EM charge dependence is constant for the 
Cg term and linear for the evCg terms, as expected from 
the smallness of the fine structure constant in QED. We 
perform uncorrelated chiral fits for the e^, esev(73-even) 
and esev(73-odd) terms separately setting /i to the chiral 
scale A^ = 1 GcV and obtain the LEC's in Tab. [Til In this 
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FIG. 3. contribution to Afps (uncorr). Lines represent 
uncorrelated fits to SU{2) PQChPT. 
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FIG. 4. Same as Fig. [3] but for ((71,93) = (+2/3,-1/3), 
sliowing the valence quark mass dependence. 
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FIG. 5. eBev(73-even) contribution to Afpg (uncorr) 
represent uncorrelated fits to SU(2) PQChPT. 
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showing the valence quark mass dependence. 



3x1 0's 
2x10 = 
1x10„ 
■9' 

-IXIO5 
-2x10' 
-3x10, 
-4x1 0's 
-5x10^ 
-6x10 



1 1 2 2 

r e^e^ (Tg-odd) term in AMpg (GeV) 
i ^ ^ X . 




- 


a (am,, am^) = (0.01, 0.01) 
A (am,, amj) = (0.02, 0.02) 
« (am,, amj) = (0.03, 0.03) 


1 




- 



1/3 



2/3 



FIG. 7 

represent 



4/3 



esev(73-odd) contribution to Mpg (uncorr) 
uncorrelated fits to SU{2) PQChPT. 



Lines 



1x10" 
5x10"' 



-5x10 
-1x10"" 



' I ' I ' r 

®s^v (Tg-odd) term in AMpg 




a('^3"^"^res(QCD)) 



FIG. 8. Same as Fig. [7] but for (gi,(/3) = ( 
showing the valence quark mass dependence. 
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fit, we choose a minimal set of data with smaller valence 
quark masses, and ignore Qq dependence in Bq because 
of smallness of and Yi. We also neglect finite vol- 
ume effects in this analysis, which could give significant 
shifts in the EM mass spectrum, since the lattice volume 
used here is small. However, we also note that the quark 
masses used here are relatively heavy. Although the sta- 
tistical error is large, the value of LEG C is consistent 
with that obtained in quenched QED 0. (The lattice 
volume and the quark masses used in the chiral fit are 
different between this work and Ref. The important 
fact, however, is that the order of magnitude is consis- 
tent between them.) The size of Yi seems to be the same 
as the other QED LEC's in ©(e^m) terms determined in 
quenched QED 0] , which means the sea EM charge effect 
is comparable to the valence one except for the Dashen 
term. 

In this study incorporating sea quark EM charges in 
2+1 flavor lattice QED+QCD, we have shown that the 
QED LEC's are accessible using the reweighting method, 
and that the sea quark LEC's are the same size as the 



TABLE 11. QED low-energy constants at the chiral scale 
= 1 GeV. 3^1 is defined as Vi = VitrQ^^, for SU{3) 
ChPT and 3^1 = yitrQ^(2) + yi'(trQ,(2))^ + Fi"g6trQ,(2) for 
SU{2) ChPT. J and K. depict J = JtrQs(2) + J'qe and 
IC — A'trQs(2) + K'qe, respectively. The quenched QED 
(qQED) values for C are quoted from Ref. whose val- 
ues are obtained from 24^ x 64 lattice and by infinite volume 
ChPT formula. The values of the LEC's Bo and Fq used in 
the chiral fit are quoted from Ref. 01 ■ 





5(7(3) ChPT 5(7(2) ChPT 
uncorr corr uncorr corr 


W'C (qQED) 
lO^C 
lO^Fi 
10^3^1 
lO^J 
10*/C 


2.2(2.0) - 18.3(1.8) 

8.4(4.3) 8.3(4.7) 20(14) 15(21) 

-5.0(3.6) -0.4(5.6) 

-3.1(2.2) -0.2(3.4) -3.0(2.2) -0.2(3.4) 
-2.6(1.6) -3.3(2.8) 
-3.1(6.9) -3.7(7.8) 



valence ones, as expected. In our analysis, the sign flip 
engineering of EM charges proved to be highly effective, 
similar to the ±e trick for the valence sector 0, Work 
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on a larger lattice, 24^ x 64, with further algorithmic im- 
provements, for example, low-mode averaging to increase 
statistics, is in progress. 
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